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Abstract 

The spin-averaged structure function of the neutron, Kf, is extracted from 
recent deuteron data, taking into account the small but significant corrections 
from nuclear effects in the deuteron. At small x, the F^/F^ ratio measured 
by the New Muon and Fermilab E665 Collaborations is interpreted to suggest 
a small amount of shadowing in deuterium, which acts to enhance K! 1 for 
x ^$ 0.1. A careful treatment of Fermi motion, binding and nucleon off-shell 
effects in the deuteron also indicates that the neutron / proton structure functio 
ratio as x — > 1 is consistent with the perturbative QCD expectation of 3/7, 
but larger than the traditional value of 1/4. 

I. INTRODUCTION 

The quark structure of the nucleon is one of the most fundamental aspects of hadron 
physics. Deep inelastic scattering (DIS) of leptons from hydrogen has yielded a wealth 
of information on the quark and gluon substructure of the proton. The absence of free 
neutron targets means, however, that it is difficult to obtain direct data on F£- As a 
result, one usually uses deuterium targets, and extracts neutron structure information from 
a knowledge of the proton structure function, and the nucleon wave function in the deuteron. 
The accuracy of the extracted neutron data naturally depends on the level of understanding 
of the nuclear physics in the deuteron, as well as on the extraction procedure itself. Both of 
these issues are carefully addressed in this paper. 



* Invited talk presented by W. Melnitchouk at the Cracow Epiphany Conference on Proton Struc- 
ture, 5-6 Jan. 1996, Cracow, Poland. To appear in Acta Phys. Pol. 
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The treatment of nuclear effects is divided into two regions: large x (x ^ 0.3) and small 
x (x ^ 0.2). In the context of the multiple scattering framework, the l&rge-x effects are 
described within the impulse approximation, in which the virtual photon interacts with only 
one nucleon in the deuteron, while the other nucleon remains spectator to the interaction. 
The impulse approximation provides a natural framework within which effects from nuclear 
binding, Fermi motion, and nucleon off-shellness can be incorporated. At small x, on the 
other hand, there are important contributions from the rescattering of the probe from both 
nucleons in the deuteron, which gives rise to the phenomenon known as nuclear shadowing. 

II. LARGE X 

Away from the small- x region (x ^ 0.3), the dominant contribution to the deuteron 
structure function can be computed from the impulse approximation. Here the total 
amplitude is factorized into 7* A and ND amplitudes, although, contrary to what is often 
assumed, this factorization does not automatically lead to a factorization of cross sections 
in the convolution model []IJ . 

A. Binding, Fermi Motion and Off-Shell Effects 

Starting from the impulse approximation, one can show that in the non-relativistic ap- 
proximation the nuclear structure function can be written in convolution form, in which 
the structure function of the nucleon is smeared with a momentum distribution, /at/d(z/)) of 
nucleons in the deuteron |T],0|: 

F 2 D (conv) (x, Q 2 ) = Jdy f N/D (y) F 2 N (j, Q 2 ) , (1) 

where F 2 N — F$ + F 2 . Equation (JJ) is correct to order iy/c) 2 (with v the nucleon velocity), 
provided we also neglect the possible p 2 dependence in the nucleon structure function. The 
distribution function In/d^v) is determined by the deuteron wave functions u,w,v t and v s 
II (corresponding to the deuteron's S, D and triplet and singlet P waves) (§]: 

f N/D (y) = ^ y f p2max d P 2 ^ (n 2 ( P 2 ) + w 2 ( P 2 ) + , 2 ( P 2 ) + , 2 ( P 2 )) (2) 

4 J-00 Po 

where p is the interacting nucleon's four-momentum, with a maximum squared value p 2 max = 
yM 2 D -yM 2 /(l-y),tmdE p = yJM 1 + p 2 . 

As explained in Refs. HH, explicit corrections to Eq.fJI]), which cannot be written in 
convolution form, arise when the bound nucleons' off-mass-shell structure is taken into ac- 
count: 

F 2 D (x, Q 2 ) = F 2 (conv) (x, Q 2 ) + 5^F 2 D (x, Q 2 ). (3) 

The correction 8^ oS ^F 2 D receives contributions from the off-shell components in the deuteron 
wave function, as well as from the off-mass-shell dependence of the bound nucleon structure 
function [|J (i.e. in the p 2 —>■ M 2 limit, in which the P-state wave functions also vanish, one 
has 5(° ff )P 2 D -> 0). 
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FIG. 1. 



/F£ ratio function of x for the off-shell model of Refs. 



(solid) and the 



on-shell model of Ref. M (dotted). 



In Refs. HH the structure function was modeled in terms of relativistic quark- 
nucleon vertex functions, which were parametrized by comparing with available data for 
the parton distribution functions. The off-shell extrapolation of the j*N interaction was 
modeled assuming no additional dynamical p 2 dependence in the quark-nucleon vertices. 
This enabled an estimate of the correction 5^ oS 'F^ to be made, which was found to be quite 
small, of the order ~ 1 — 2% for x ^ 0.9. The result of the fully off-shell calculation from Ref. 
|| is shown in Fig.l (solid curve), where the ratio of the total deuteron to nucleon structure 
functions {F^* / F% ) is plotted. Shown also is the result of an on-mass-shell calculation from 
Ref. (dotted curve), which has been used in many previous analyses of the deuteron data 
The most striking difference between the curves is the fact that the on-shell ratio has 
a very much smaller trough at x ~ 0.3, and rises faster above unity (at x ~ 0.5) than the 
off-shell curve, which has a deeper trough, at x ~ 0.6 — 0.7, and rises above unity somewhat 
later (at x ~ 0.8). 

The behavior of the off-shell curve in Fig.l is qualitatively similar to that found by 
Uchiyama and Saito ||, Kaptari and Umnikov ||10|| , and Braun and Tokarev [ITJ, who also 
used off-mass-shell kinematics, but did not include the (small) non-convolution correction 
term S^'F^ ■ The on-shell calculation 0], on the other hand, was performed in the infinite 
momentum frame where the nucleons are on their mass shells and the physical structure 
functions can be used in Eq. (|TJ) . One problem with this approach is that the deuteron 
wave functions in the infinite momentum frame are not explicitly known. In practice one 
usually makes use of the ordinary non-relativistic S- and D-state deuteron wave functions 
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calculated in the deuteron rest frame, a procedure which is analogous to including only 
Fermi motion effects in the deuteron. In addition, one knows that the effect of binding in 
the infinite momentum frame shows up in the presence of additional Fock components (e.g. 
iViV-meson(s) ) in the nuclear wave function, which have not yet been computed. 

Clearly, a smaller D/N ratio at large x, as in the off-shell calculation, implies a larger 
neutron structure function in this region. To estimate the size of the effect on the n/p ratio 
requires one to "deconvolute" Eq.(j^) in order to extract F£. 



B. Extraction of F£ 

To study nuclear effects on the neutron structure function arising from different models 
of the deuteron, one must eliminate any effects that may arise from the extraction method 
itself. We therefore use exactly the same extraction procedure as used in previous SLAC 
and EMC M data analyses, namely the smearing (or deconvolution) method discussed 



by Bodek et at. ||12|| . This method involves the direct use of the proton and deuteron data, 
without making any assumption for F% itself. For completeness let us briefly outline the 
main ingredients in this method. (For alternative methods of unfolding the neutron structure 
function see for example Refs. jO,|TJ|.) 

Firstly, one subtracts from the deuteron data, Ff, the additive, off-shell corrections, 
5( oS )F®, to give the convolution part, F® (conv ). Then one smears the proton data, F$, with 
the nucleon momentum distribution function f^m{y) in Eq.([l|) to give F% = F% / S p . The 
smeared neutron structure function, F^, is then obtained from 

pn = p D (conv) _ f p ^ 

Since the smeared neutron structure function is defined as F£ = F£ / S ni we can invert this 
to obtain the structure function of a free neutron, 

F? = S n (F 2 D (conv) - F*/S p ) . (5) 

The proton smearing factor, S p , can be computed at each x from the function fxm{y), 
and a parametrization of the data (for example, the recent fit in Ref. [15] to the combined 



SLAC, BCDMS and NMC data). The neutron F£ structure function is then derived from 
Eq.(||) taking as a first guess S n = S p . These values of F% are then smeared by the function 
In/d^v), and the results used to obtain a better estimate for S n . The new value for S n is 
then used in Eq.(0) to obtain an improved estimate for F£, and the procedure repeated until 
convergence is achieved. 
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FIG. 2. F 2 /F 2 ratio as a function of x, for the off-shell model (solid), off-shell model without 
the convolution-breaking term (dashed), and the on-shell model (dotted). On the right-hand axis 
are marked the x — ► 1 limits of the SU(6) symmetric model (2/3), and the predictions of the models 
of Refs. fHH (1/4) and ©P (3/7). 

The results of this procedure for F 2 /F 2 are presented in Fig.2, for both the off-shell 
calculation (solid) and the on-shell model (dotted). The increase in the off-shell ratio at 
large x is a direct consequence of the deeper trough in the F 2 /F 2 ratio in Fig.l. To 
illustrate the role of the non-convolution correction, S ( -°^F 2 D ', we have also performed the 
analysis setting this term to zero, and approximating F 2 by F 2 ^ c0 (x). The effect of this 
correction (dashed curve in Fig.2) appears minimal. One can therefore attribute most of the 
difference between the off- and on-shell results to kinematics, since both calculations involve 
very similar deuteron wave functions. 

Before discussing the implications of these results, let us briefly outline the connection 
between structure functions at large x and the breaking of SU(6) spin-flavor symmetry. 



C. SU(6) Symmetry Breaking 

The large-x region, being valence quark dominated, is where SU(6) symmetry breaking 
effects in valence quark distributions should be most prominent. The precise mechanism for 
the breaking of the spin-flavor SU(6) symmetry is a basic question in hadronic physics. In 
a world of exact SU(6) symmetry, the wave function of a proton, polarized say in the +z 
direction, would be simply [115): 



V T = 77=" T {ud) s =o + -j=u T (ud) s =i ~ \ u I ( ud )s=i 
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- id t (uu) s =i - i (uu) s =i, (6) 

where the subscript 5 denotes the total spin of the two-quark component. In this limit, 
apart from charge and flavor quantum numbers, the u and d quarks in the proton would be 
identical. The nucleon and A isobar would, for example, be degenerate in mass. In deep- 
inelastic scattering, exact SU(6) symmetry would be manifested in equivalent shapes for the 
valence quark distributions of the proton, which would be related simply by uv(x) = 2dy(x) 
for all x. For the neutron to proton structure function ratio this would imply: 



F 



[SU(6) symmetry]. (7) 



In nature spin- flavor SU(6) symmetry is, of course, broken. The nucleon and A masses 
are split by some 300 MeV. Furthermore, with respect to DIS, it is known that the d quark 
distribution is softer than the u quark distribution, with the neutron/proton ratio deviating 
at large x from the SU(6) expectation. The correlation between the mass splitting in the 56 



baryons and the large-x behavior of Fg/F^ was observed some time ago by Close |T7[ and 



Carlitz ||18|| . Based on phenomenological [fL7] and Regge |T8[ arguments, the breaking of the 
symmetry in Eq.(|6]) was argued to arise from a suppression of the "diquark" configurations 
having 5 = 1 relative to the 5 = configuration, namely 

(qq)s=o > (qq)s=i, x->l. (8) 

Such a suppression is in fact quite natural if one observes that whatever mechanism leads 
to the observed N — A splitting (e.g. color-magnetic force, instanton-induced interaction, 
pion exchange), it necessarily acts to produce a mass splitting between the two possible spin 
states of the two quarks, (qq)s, which act as spectators to the hard collision, with the 5 = 1 



state heavier than the 5 = state by some 200 MeV [19]. From Eq.(^), a dominant scalar 
valence diquark component of the proton suggests that in the x — > 1 limit F$ is essentially 
given by a single quark distribution (i.e. the it), in which case: 

F n 1 d 

-^p — > -, > [5 = dominance]. (9) 

This expectation has, in fact, been built into many phenomenological fits to the parton 
distribution data. 

An alternative suggestion, based on perturbative QCD, was originally formulated by 
Farrar and Jackson ||20|| . There it was argued that the exchange of longitudinal gluons, 
which are the only type permitted when the spins of the two quarks in (qq)s are aligned, 
would introduce a factor (1 — x) 1 / 2 into the Compton amplitude — in comparison with the 
exchange of a transverse gluon between quarks with spins anti-aligned. In this approach the 
relevant component of the proton valence wave function at large x is that associated with 
states in which the total "diquark" spin projection, S z , is zero: 

(g?k=o > (w)s,=i, (io) 

Consequently, scattering from a quark polarized in the opposite direction to the proton 
polarization is suppressed by a factor (1 — x) relative to the helicity-aligned configuration. 
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A similar result is also obtained in the treatment of Brodsky et al. (based on 
counting- rules) , where the large- x behavior of the parton distribution for a quark polar- 
ized parallel (AS Z = 1) or antiparallel (AS Z = 0) to the proton helicity is given by: 
q^(x) = (1 — x) 2n_1+A5z , where n is the minimum number of non- interacting quarks 
(equal to 2 for the valence quark distributions). In the x — > 1 limit one therefore predicts: 

^§r -> - -> \ [S z = dominance]. (11) 

Note that the d/u ratio does not vanish in this model. Clearly, if one is to understand the dy- 
namics of the nucleon's quark distributions at large x, it is imperative that the consequences 
of these models be tested experimentally. 

The reanalyzed SLAC ||7|,p2| data points themselves are plotted in Fig. 3, at an average 
value of Q 2 ~ 12 GeV 2 . The very small error bars are testimony to the quality of the SLAC p 
and D data. The data represented by the open circles have been extracted with the on-shell 
deuteron model of Ref. 0, while the filled circles were obtained using the off-shell model of 
Refs. Most importantly, the F^jF^ points obtained with the off-shell method appear 



to approach a value broadly consistent with the Farrar- Jackson [£(J and Brodsky et al. pi 
prediction of 3/7, whereas the data previously analyzed in terms of the on-shell formalism 
produced a ratio that tended to the lower value of 1/4. 
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FIG. 3. Deconvoluted F 2 n /F 2 p ratio extracted from the SLAC p and D data g^, at an average 
value of Q 2 ~ 12 GeV 2 , assuming no off-shell effects (open circles), and including off-shell effects 
(full circles). 

The d/u ratio, shown in Fig. 4, is obtained by inverting Fg/F^ in the valence quark 
dominated region. The points extracted using the off-shell formalism (solid circles) are 
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again significantly above those obtained previously with the aid of the on-shell prescription. 
In particular, they indicate that the d/u ratio may actually approach a finite value in the 
x — > 1 limit, contrary to the expectation of the model of Refs. |OJT3], in which d/u tends 
to zero. Although it is a priori not clear at which scale the model predictions |l7|JT8| . [20| . [2l[| 
should be valid, for the values of Q 2 corresponding to the analyzed data the effects of Q 2 
evolution are minimal. 
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FIG. 4. Extracted d/u ratio, using the off-shell deuteron calculation (full circles) and using 
on-shell kinematics (open circles). Also shown for comparison is the ratio extracted from neutrino 
measurements by the CDHS collaboration [p3|]. 

Naturally it would be preferable to extract F% at large x without having to deal with 
uncertainties in the nuclear effects. In principle this could be achieved by using neutrino 
and antineutrino beams to measure the u and d distributions in the proton separately, and 
reconstructing F^ from these. Unfortunately, as seen in Fig. 4, the neutrino data from the 
CDHS collaboration J23| do not extend out to very large x (x ^ 0.6), and at present cannot 
discriminate between the different methods of analyzing the electron-deuteron data. 



The results of our off-shell model are qualitatively similar [22J to those obtained using the 

There the EMC effect 



nuclear density method suggested by Frankfurt and Strikman |24 



in deuterium was assumed to scale with that in heavier nuclei according to the ratio of the 
respective nuclear densities, so that the ratio F 2 D /F^ V in the trough region was depleted by 
about 4%, similar to that in Fig.l (solid curve). This would give an F^/F^ ratio broadly 
consistent with 3/7. 

We should also point out similar consequences for the spin-dependent neutron structure 
function g™, where the models of Refs. flr?| , |IB] and Refs. also give different predictions 

for g™/gi as x — > 1, namely 1/4 and 3/7, respectively. Quite interestingly, while the ratio of 



S 



polarized to unpolarized u quark distribution is predicted to be the same in the two models, 

Am 



it 



1 [S = or S z = dominance], (12) 



the results for the <i-quark distribution ratio differ even in sign: 

— > — — [S = dominance], (13a) 

— > 1 [S z = dominance]. (13b) 

To extract information on the polarized parton densities at large x that is capable of dis- 
criminating between these predictions, the same care will need to be taken when subtracting 
the nuclear effects from g^ and g\ He . In particular, the results of Refs. PJ25|| indicate that 



while the simple prescription [^B|] of subtracting the g\ structure function from the D data, 
modified only by the deuteron D-state probability, is surprisingly good for x ^ 0.6, it is 
completely inadequate for x ^ 0.7. 

Having seen that how one handles nuclear corrections can critically affect the deciphering 
of the physical implications of the extracted structure function at large x, we now examine 
the consequences of nuclear effects in deuterium in the small- a; region. 



III. SMALL X 

At small values of x the impulse approximation should eventually break down. Indeed, 
one finds that coherent multiple scattering effects become very important when the charac- 
teristic time scale 1/Mx of the DIS process becomes larger than the typical average distance 
between bound nucleons in the nucleus, which occurs typically for x ^ 0.1. These effects are 
seen, for example, in the low-x depletion of the nuclear EMC ratio, F 2 / F 2 D . Any shadowing 
in the deuteron itself should therefore produce a depletion in the F 2 / F 2 ratio at small x 

2%m\. 



A. Nuclear Shadowing 

The rescattering of the virtual photon from several nucleons in a nucleus is usually 
described within the non-relativistic Glauber multiple scattering formalism. Relativistic 
corrections will amount to a few percent out of a shadowing correction to F 2 that will itself 
be a few percent in total, and hence can be neglected. For the deuteron the only contribution 
in the Glauber series comes from the double scattering process. 

At small x, nuclear binding and Fermi motion corrections can be neglected, and the total 
deuteron structure function written as: 

F 2 D (x, Q 2 ) » Fi(x, Q 2 ) + F?(x, Q 2 ) + 6^F 2 D (x, Q 2 ). (14) 

In modeling the shadowing correction, 5 < - shad ^ F 2 D , our approach is to take a two-phase model, 
similar to that of Kwiecinski and Badelek p7|-p9| . At high virtuality the interaction of the 
virtual photon with the nucleus is parametrized in terms of diffractive scattering through 
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the double and triple Pomeron, as well as scattering from exchanged mesons in the deuteron. 
On the other hand, at low virtuality it is most natural to apply a vector meson dominance 
(VMD) model, in which the virtual photon interacts with the nucleons via its hadronic 
structure, namely the p°, u and mesons. The latter contribution vanishes at sufficiently 
high Q 2 , but for Q 2 < 1 GeV 2 it is in fact responsible for the majority of the Q 2 variation. 

For the diffractive component, Pomeron (IP) exchange between the projectile and two or 
more constituent nucleons models the interaction of partons from different nucleons within 
the deuteron. Assuming factorization of the diffractive cross section, the shadowing correc- 
tion (per nucleon) to the deuteron structure function from iP-exchange is written as 
a convolution of the Pomeron structure function, -F 2 /p , with a distribution function ("flux 
factor"), fip/D, describing the number density of exchanged Pomerons: 

6^F° (x, Q 2 ) = f dy f IP/D {y) F?(x IP , Q 2 ), (15) 



where, within the non-relativistic approximation f27H3l 



fi P/D (y) = -^\j dV S D (k 2 ). (16) 



Here y = x(l + M x /Q 2 ) the light -cone momentum fraction carried by the Pomeron (Mx is 
the mass of the diffractive hadronic debris), and xrp = x/y is the momentum fraction of the 
Pomeron carried by the struck quark in the Pomeron. The deuteron form factor, Sn(h 2 ), is 
given in terms of the coordinate space wave functions: 



Sd(V) = J q dr[u 2 (r) + w 2 (r))j (\k\r), (17) 

where jo is the spherical Bessel function. Within experimental errors, the factorization hy- 
pothesis, as well as the y dependence of fip/A{y) H3~Hl> are consistent with the recent HERA 
data |33| obtained from observations of large rapidity gap events in diffractive ep scattering. 
These data also confirm previous findings that the Pomeron structure function contains both 
a hard and a soft component: F 2 p (a;/p, Q 2 ) = F 2 /P (xip, Q 2 ) + F2 P ^° h \xip, Q 2 ). The hard 
component of _F 2 /P is generated from an explicit qq component of the Pomeron, and has an 
xip dependence given by x/p(l — £/p) in agreement with the recent diffractive data 
33fl . The soft part, which is driven at small Xj P by the triple- Pomeron interaction ||27|| , has 



a sea quark-like xjp dependence, with normalization fixed by the triple-Pomeron coupling 
constant. 

The dependence of F 2 P on Q 2 at large Q 2 arises from radiative corrections to the par- 
ton distributions in the Pomeron [28], which leads to a weak, logarithmic, Q 2 dependence 
for the shadowing correction 8( IpS) F®. The low-Q 2 extrapolation of the qq component is 
parametrized by applying a factor Q 2 /(Q 2 + Q%), where Ql f& 0.485 GeV 2 |35| may be in- 



terpreted as the inverse size of partons inside the virtual photon. For the nucleon sea quark 
densities relevant for j?^ p ^ oit ^ we use the recent parametrization from Ref. ||35|| , which in- 
cludes a low-Q 2 limit consistent with the real photon data, in which case the total Pomeron 
contribution 5^ IP) F^ -> as Q 2 -> 0. 

To adequately describe shadowing for small Q 2 requires one to use a higher-twist mech- 
anism, such as vector meson dominance. VMD is empirically based on the observation that 
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some aspects of the interaction of photons with hadronic systems resemble purely hadronic 
interactions. In terms of QCD this is understood in terms of a coupling of the photon to 
a correlated qq pair of low invariant mass, which may be approximated as a virtual vector 
meson. One can then estimate the amount of shadowing in terms of the multiple scattering 
of the vector meson using Glauber theory. The corresponding correction (per nucleon) to 
the nuclear structure function is: 



where 

Sa VD 



a 



VN 



J d 2 k T S D (k 2 ) (19) 



is the shadowing correction to the vector meson — nucleus cross section, fy is the photon — 
vector meson coupling strength, and My is the vector meson mass. 

In practice, only the lowest mass vector mesons (V = p°,u,<p) are important at low 
Q 2 . For Q 2 — > and fixed x, S^F® disappears due to the vanishing of the total ■ 
Furthermore, since this is a higher twist effect, shadowing in the VMD model dies off quite 
rapidly between Q 2 ~ 1 and 10 GeV 2 , so that for Q 2 > 10 GeV 2 it is almost negligible - 
leaving only the diffract ive term, S^'F^ ■ (Note that at fixed is, for decreasing Q 2 the ratio 
F® '/F$ approaches the photoproduction limit.) 

Another potential source of shadowing arises from the exchange of mesons between the 
nucleons and the probe. It has previously been suggested [TJ(J that this leads to some anti- 
shadowing corrections to F®(x). The total contribution to the deuteron structure function 
from meson exchange is: 



, K x, Q 2 ) = W dy f M/D (y) F 2 M (x/y, Q 2 ), (20) 

» r J X 



M 



where M = ir,p,uj,a. The virtual meson structure function, F 2 M , one approximates by 
the (real) pion structure function, data for which exist from Drell-Yan production. The 
exchange- meson distribution functions fuiu) are obtained from the non-relativistic reduc- 



tion of the nucleon — meson interaction given in Refs. |30|j36[| . In practice pion exchange is 
the dominant process, and this gives a positive contribution to 8^ M ^ > F 2 D {x ) Q 2 ). The exchange 
of the fictitious a meson (which represents correlated 2ir exchange) also gives rise to anti- 
shadowing for small x. Vector mesons (p, uj) exchange cancels some of this antishadowing, 
however the magnitude of these contributions is smaller. In fact, for soft meson-nucleon 
vertices (Am ~ 1.3 GeV) all contributions other than that of the pion are totally negligible. 



B. Neutron Structure Function at Small x 

For x ^ 0.1 the magnitude of the (negative) Pomeron/VMD shadowing is larger than the 
(positive) meson-exchange contribution, so that the total #( shad ).F 2 D is negative. For larger 
x (~ 0.1 — 0.2) there is a very small amount of antishadowing, which is due mainly to the 
VMD contribution, and also to the pion-exchange contribution. For the NMC kinematics 
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(x > 0.004, Q 2 = 4 GeV 2 ) f37fl , the overall effect on the shape of the neutron structure 
function is a 1 — 2% increase in _pny_p"> bound f or x < 0.01, where ^r™> bound = pD _ pP m 
The presence of shadowing in the deuteron would be confirmed through observation of a 
deviation from unity in the F^/F^ structure function ratio in the kinematic region where 
Regge theory is expected to be valid. Although the exact value of x below which the proton 
and (free) neutron structure functions become equivalent is not known, it is expected that 
at low enough x, if -> F%, in which case F 2 D /F| -> 1 + 5^ & ^F°/F%. While for the 
lowest NMC data point it may be debatable whether the Regge region is reached, the E665 
Collaboration has taken data to very low x, x ~ 10~ 5 , which should be much nearer the 
onset of Regge behavior. 
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FIG. 5. x dependence of the D/p structure function ratio, compared with the low-x E665 
data |38|] and NMC data [37] at larger x. The dashed curve is the result without any shadowing 
correction. 



In Fig. 5 we show the low x E665 Collaboration data ||38|| , as well as the earlier NMC 
data at larger x p7| (note that the E665 data are not taken at fixed Q 2 ). The calculated 
ratio with a small shadowing correction is shown by the solid curve, while the result for the 
case of no shadowing is indicated by the dashed curve. Assuming F% = F£ in the lowest 
x-bins, the data clearly favor the shadowing scenario. 



C. Are 1% Effects Worth Worrying About? Gottfried Sum Rule 



The main interest in the NMC measurment of the neutron structure function at low x 
was to test accurately the Gottfried sum rule, 
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Sc= t dx m^im (21a) 

Jo X 

= -- f dx (d(x) - u(x)) (21b) 
3 Jo 

which, in the naive quark model where d = u, is S G m = 1/3. Ignoring nuclear effects, the 
experimental value obtained by the NMC was S^ p = 0.258 ± 0.017, indicating a violation 
of SU(2) flavor symmetry in the proton's sea. However, because the structure function 
difference in Sg is weighted by a factor 1/x, any small differences between F% and p^ hound 
could be amplified for x — ► 0. 

Including the above shadowing correction, the overall effect on the experimental value 
for S G , 

Sg = S^+ / dx- (22) 

JO X 

is a reduction of between -0.010 and -0.026, or about 4 and 10% of the measured value 
without putting in the shadowing correction |5U] . Therefore a value that reflects the "true" 



Gottfried sum should be around Sg ~ 0.2. This is some 30% reduction from the naive quark 
model prediction, S G m . Hence the tiny shadowing effect in the deuteron results in a more 
significant violation of flavor SU(2) symmetry in the proton sea. 



IV. SUMMARY 

In summary, we have reanalyzed the latest proton and deuteron structure function data 
at small and large x, in order to obtain more reliable information on the structure of the 
neutron in the x — > and x — > 1 limits. 

At small x, we have estimated the nuclear shadowing arising from the double scattering 
of the virtual photon from both nucleons in the deuteron. To cover the entire range of Q 2 
accessible to current experiments, the 7*iV interaction is described in terms of the VMD 
model, which models the correlated qq pair excitations of the virtual photon, together with 
Pomeron exchange for the uncorrelated qq pair interactions with the nucleon. Such a hybrid 
model is particularly necessary if one is to reproduce the observed Q 2 dependence of the EMC 
ratio at small x |39|j . In addition, we have also included contributions from the exchange of 
mesons, which effectively cancels as much as half of the shadowing from the VMD/Pomeron- 
exchange mechanisms alone. The net effect is a ^ 1% reduction of F 1 ? for x ~ 0.004, or 
equivalently a ^ 2% increase in the neutron structure function over the uncorrected F£. 



This is consistent with the recent measurement by the Fermilab E665 Collaboration |38| of 
the ratio F® jF% down to x ~ 10 -5 , which deviates from unity by about 4%, and suggests 
the presence of shadowing in the deuteron. Although the absolute values of the shadowing 
corrections to FJ? are small, because they are concentrated at small x, the effect on the 
Gottfried sum S G is a further reduction of up to ~ 10% over the value measured by NMC 

Including all of the currently known nuclear effects in the deuteron at large x, namely 
Fermi motion, binding, and nucleon off-mass-shell effects, we find that the total EMC effect 
is ~ 2% larger than in previous calculations based on on-mass-shell kinematics, from which 
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binding effects were omitted. The larger deviation from unity for 0.5 ^ x ^ 0.8 in /FJ? 
translates into an increase in the ratio Ff /Ff. Our results indicate that as x — > 1 the limiting 
value of F^/Ff is above the previously accepted result of 1/4, and broadly consistent with 
the perturbative QCD expectation of 3/7. This also implies that the d/u ratio approaches 
a non-zero value of 1/5 as x — > 1 . 

Finally, for definitive tests of the nuclear effects in the deuteron one would like model- 
independent information on the neutron at both low and high x. In principle, this can 
be achieved with high-precision data from neutrino-proton experiments, from which indi- 
vidual flavor distributions can be determined, and the neutron structure function inferred 
from charge symmetry. Unfortunately, both the statistics on the neutrino data and the 
coverage in x will not allow this in the near future. It has therefore been suggested that 
one might perform a series of semi-inclusive experiments on deuteron targets, measuring 
in coincidence both the scattered lepton and recoiling proton or neutron. This would help 
constrain the deuteron wave function over a large range of kinematics, and hence enable the 
^/-dependence of the nucleon momentum distribution functions to be mapped out directly. 
Such experiments are already planned for CEBAF and HERMES |3(J, and should provide 
critical information on the size and importance of relativistic and other short- distance nu- 
clear phenomena in the deuteron. 
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